Random walk in a two-dimensional self-affine random potential : 
properties of the anomalous diffusion phase at small external force 
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We study the dynamical response to an external force F for a particle performing a random walk 
in a two-dimensional quenched random potential of Hurst exponent H — 1/2. We present numerical 
results on the statistics of first-passage times that satisfy closed backward master equations. We 
find that there exists a zero- velocity phase in a finite region of the external force < F < F c , where 
the dynamics follows the anomalous diffusion law x(i] ~ £(F) t^ F \ The anomalous exponent 
< fJ,{F) < 1 and the correlation length £(F) vary continuously with F. In the limit of vanishing 
£\J ' force F — >• 0, we measure the following power-laws : the anomalous exponent vanishes as fi(F) oc F a 

with a ~ 0.6 (instead of o = 1 in dimension d — 1), and the correlation length diverges as £(F) tx F~" 
with v ~ 1.29 (instead of v — 2 in dimension d = 1). Our main conclusion is thus that the dynamics 
renormalizes onto an effective directed trap model, where the traps are characterized by a typical 
length £(F) along the direction of the force, and by a typical barrier l/fx(F). The fact that these 
^ , traps are 'smaller' in linear size and in depth than in dimension d = 1, means that the particle uses 

m 



bo- 
3 



the transverse direction to find lower barriers 



Ph ! I. INTRODUCTION 

In quantum dynamics, the discovery of Anderson localization fifty years ago [l[ has shown that disorder can 
completely suppress transport over large distances. For classical random walks, the presence of quenched disorder can 
. also induce completely new transport behaviors with respect to pure systems in some regions of parameters (see the 
reviews 043 and references therein). In the present paper, we study the dynamical response to an external force F 
Q f° r a particle performing a random walk in a two-dimensional quenched random potential U(r). The dynamics of the 
' position r{t) of the particle is governed by the Langevin equation 

Ch ' dr 

O: -=F-VU(f(t)) + ff(t) (1) 



dt 

where ff(t) is a white noise 



*> ■ < Vi{t)Vj(t') >= 2T8(t - t')6 id (2) 

' and where the quenched random potential U(r) is self-affine with some Hurst exponent H 



[U{r)-U{r')] 2 ^ \r-r'\ 2H (3) 

I r — r ' I — > oo 

The random-force Sinai model that corresponds to dimension d = 1 and Hurst exponent H = 1/2, has been much 
studied by various exact methods (see for instance the review [ij and references therein) with the following conclusions: 

(i) for the case F = 0, the dynamics is logarithmically-slow \x(t) — x(0)\ ~ (\nt) 2 ; 

(ii) for the case F ^ 0, there exists various dynamical phase transitions as a function of F = \F\ [^-Q. The 
first important transition concerns the velocity V defined as the ratio between the thermally averaged displacement 

' < x {t) > an d the time t in the limit t +oo 



< x(t) > ~ V(F)t + ... for F > F c 

< x(t) > ~ t^ F) X + ... for F < F c (4) 

For F > F c , the velocity V(F) is finite and self-averaging with V(F — > F c ) — 0. For F < F c , the velocity vanishes 
V = 0, and the leading term of the displacement is governed by a continuously varying exponent < fJ.(F) < 1 with 
fx(F — > 0) = and fj,(F — > F c ) — > 1. In addition, the reseated variable X =< x(t) > /t^ F ^ remains distributed over 
the samples. In this anomalous diffusion phase, the width of the diffusion front in a given sample in also of order 

(<x 2 (t) > - < x(t) > 2 ) 1/2 ~f( F) for F<F C (5) 

More detailed statements about the diffusion fronts can be found in the review [ij . From the point of view of strong 
disorder renormalization fol— fl~lT| . the logarithmic behavior x(t) ~ (lnt) 2 for F = corresponds to an "infinite disorder 
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fixed point", whereas the anomalous diffusion phase x(t) ~ t^ F ^ corresponds to a "strong disorder fixed point" in 
the limit of small external force F — >• where fi(F) — > 0. 

In higher dimension d > 1, exact results are not available anymore, but the logarithmic scaling at zero external 
force is still expected from scaling arguments on barriers 0, [H| 

|r(t) -f(0)| - (lnt) 1 /^ for ^ = (6) 



We have recently checked numerically this behavior by various methods in dimension d — 2 [13H15I ]. In the review 
[H, the response to an external force F in dimension d > 1 was left as an open question. We are not aware of works 
on this question since. The aim of the present paper is thus to study numerically the response to an external force 
in dimension d — 2, and to determine in particular whether there exists a zero- velocity phase as in dimension d = 1. 
We will focus on observables that satisfy closed backwards equations and their corresponding exact renormalization 
rules . The fact that first- passage times (or similar observables like escape probabilities through a given boundary) 
satisfy 'backwards master equation' is of course very well-known and can be found in most textbooks on stochastic 
processes (see for instance (l6l - [l9l |). In the field of disordered systems, the backward Fokker-Planck equation has 
been very much used to characterize the dynamics of a single particle in a random medium (see for instance (2p| - [25j ). 
but to the best of our knowledge, this approach has not much been used in higher dimension, nor for many-body 
problems. In our recent work [15|, we have presented numerical results for a mean-field spin model, and the present 
two-dimensional random walk model without external force, i.e. for the case F = in Eq. [TJ The present paper is 
devoted to the case of a non-zero external force F > 0. 

The paper is organized as follows. In section [TTl we define more precisely the model used in our numerical study. 
Section ITO1 concerns escape probabilities, whereas section|lV]is devoted to first-passage times. In section IVl we discuss 
the interpretation of our results in terms of a renormalized directed trap model, and we discuss the similarities and 
differences with the same model in dimension d = 1. Our conclusions are summarized in section [VTl 



We consider the following master equation that describes the time evolution of the probability P t (r) to be at position 
r at time t 



II. MASTER EQUATION DEFINING THE STOCHASTIC DYNAMICS 

We consider a two-dimensional square lattice of size L x L. To generate a quenched random self-affine potential 
U{r) of Hurst exponent H = 1/2 (Eq[3]), we have used the so-called Weierstrass-Mandelbrot function method (see 
131 and references therein for more details). In the presence of an external force F in the ^-direction, it is useful to 
introduce the tilted potential Up defined as 

U F (x,y) = U(x,y)-Fx (7) 

,tion ' 

= J2 P t (r ') W{r' ^ f) - P t (r) W out (r) (8) 

r ' 

where W (r ' — > r) represents the transition rate per unit time from position r f to r, and 

W out {r) = Y,W{r^r') (9) 

r' 

represents the total exit rate out of position f. We have chosen the Metropolis dynamics at temperature T — 1 defined 
by the transition rates 

W(?-H f ')=S < ? t ?. > min(l,e-( UF ^')-uH?))/T\ (10) 

in terms of the tilted potential of Eq. [7] The first factor S^.p > > means that the two positions are neighbors on 
the two-dimensional lattice. Various numerical results concerning the case F = can be found in our recent works 
[l3l - [l5l |. In the following, we consider only the case F > 0. 



III. STATISTICS OF THE ESCAPE PROBABILITY AGAINST THE EXTERNAL FORCE 

A. Backward master Equation for the escape probabilities 

In this section, we consider the dynamics defined by the Master Equation of Eq. Hon a square lattice of size L x L, 
with absorbing boundaries along the lines x = 1 and along x — L, and reflecting boundaries along the lines y = 1 and 
along y = L. 
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FIG. 1: We consider the lattice continuous random walk in the random potential containing an external force F in the x- 
direction (Eq. [7|). The escape probability E(x,y), defined as the probability to reach first the boundary at x = 1 before the 
boundary at x = L, when starting at position (x,y), satisfies the backward master equation of Eq. [TTJwith the boundary 
conditions of Eq. If 21 Below, we give numerical results on the typical behavior of E cen ter = E(x = L/2, y = L/2) as a function 
of the size L and of the external force F. 

The probability E(r) to reach first the absorbing boundary at x = 1 (before the other absorbing boundary at 
x = L), when starting at position r, satisfies the 'backward master equation' [IH - HH 



Y^W(r -^f')E(f')-E{r)W out {f)=0 



(11) 



for 2 < x < L — 1 and 1 < y < L, with the boundary conditions 

E(x = l,y) = 1 for 1 < y < L 

E(x = L,y) =0 for 1 < y < L (12) 

as illustrated on Fig. [T] 

As explained in 15], if one eliminates iteratively the sites, the renormalized equations for the remaining sites keep 
the same form as Eq. [Tl] with renormalized rates 

]T W R (r ^t')E (f ') - (f) E (f) = (13) 



Upon the elimination of the site tq via 



the renormalization rules read 

w R,new f ') = W R + 

In the following, we focus on the escape probability 



W R (r -> f ) W R (f -H- f ') 
W R t (f ) 
W R (f -> f ) W R (r -j. r) 
W R t (fo) 



(14) 



(15) 
(16) 



Ecenter — -£/(?" center*) 

when starting at the center r ce „t er = (x = L/2,y = L/2) of the square, in order to characterize the asymmetry 
introduced by the external force F. We may thus eliminate all sites in the region 2 < x < L — 1 except the single site 
r C enter = (x = L/2,y = L/2). Using Eq. fTJl we finally obtain 

1 



-^center ~ ^(j'center') 



W R t (reenter) 4: 



W R (r center ~> r S ) 



(17) 



in terms of the set of points fg = (a; = 1, 1 < y < L) on the line x — 1 where E(rs) — 1 
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B. Numerical results on escape probabilities 




FIG. 2: (Color on line) Typical behavior of the escape probability E cen t er defined in Eq. [16] : (a) For each value of the force 
0.1 < F < 1, the typical value -E'enter decays exponentially with the size L (Eq. I19[) , (b) Inverse of the typical correlation 
length £,typ(F) (Eq. I19[) as a function of the force F : we obtain the linear behavior £t yP (F) = 1/F. 

For various values of the force 0.1 < F < 1, we show on Fig. [2] (a) the exponential decay with the size L of the 
typical value 



ElZ ter ee e* E ~ (18) 
This defines a typical correlation length ^typ(F) via 

ln£^ ter = ln£ center ~ + ... (19) 

L^oo £typ(-t< ) 

that will diverge in the limit of vanishing force 

^ F K%oF^ (20) 

As shown on Fig. [5] (b), the typical exponent is 

v typ = 1 (21) 
We find actually that even far from the F — > limit, the typical correlation length is exactly given by 

Zty P (F) = i (22) 

This means that the correlation length seen in the escape probability E center directly reflects the linear asymmetry 
introduced by the force in the tilted random potential of Eq. [7] So this length scale has exactly the same expression 
as in the one-dimensional case 0, [ll[ and in equivalent quantum spin chains [TP . l26j . 

However, it is well known that in one dimension there exists two different correlations lengths that diverge with 
different exponents as F —> [ll], [H[ . In the present section, we have found that the typical correlation length is 
governed by the same exponent v tyv = 1 as in d = 1. In the following section, we will study the other correlation 
length £(F) that appears in the statistics of first-passage times. 
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IV. STATISTICS OF FIRST-PASSAGE TIMES 
A. Backward master equation for the first-passage times 

Suppose the dynamics starts at time t = in position r, and one is interested in the random time t where the 
dynamics will reach for the first time any point on the line x — L. As is well known (see for instance the textbooks 
|16l - [l9j ). the mean first- passage time r(f) =< t > (where the notation < . > represents the average with respect to 
the dynamical trajectories) satisfies the following 'backward master equation' for all points outside the line x = L 



Y,W(r f ') T (f ') - W out (f) r (f) = -1 



(23) 



whereas points all on the line x = L satisfy the boundary condition 

t(x = L,y) = for 1 < y < L 

as illustrated on Fig. [3] 

y ^ 



(24) 
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FIG. 3: The mean first-passage time r(x,y) on the boundary x = L when starting at position (x,y), satisfies the backward 
master equation of Eq. [23] with the boundary conditions of Eq. 1241 Below, we give numerical results on the typical behavior 
of the crossing time T crosa — t(x = 1, y = L/2) as a function of the size L and of the external force F. 



Again as explained in [15|, if one wishes to eliminate iteratively the sites, the renormalized equation will be of the 
form similar to Eq. [23] but with coefficients different from ( — 1) on the right hand-side 



E W R (r ^r')r (r ') - W R ut (r) r (f) = -K R (f 

r 1 

Upon the elimination of the site fo 



(25) 



T (ro) 



1 



W* t (fb) 



J2 wR (ro -^r")r{r")+K R {r ) 



(26) 



the renormalization rules for the rates W R are the same as in Eq. 1151 whereas the new renormalization rule for the 
inhomogeneous term of the backward master equation reads 



K 



R.new 



(r) = K R (r) + 



W R (r -> f ) 
W R t (fo) 



K R (f ) 



(27) 



Here we wish to focus on the crossing time 



Tcross =r(x = l,y = L/2) 



(28) 
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representing the mean averaged time to reach any point of the line x — L when starting in the middle (y = L/2) of 
the opposite line x = 1. Numerically, one can eliminate all points x < L — 1 except the point ri ast = (x = 1, y = L/2). 
The crossing time is then computed as 



r(x = l,y = L/2) 



K R (x = l,y = L/2) 
W D R ut (x = l,y = L/2) 



(29) 



B. Numerical results on first-passage times 
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FIG. 4: (Color on line) Typical behavior of the crossing time T cross defined in Eq. [28] : (a) For each value of the force 
O.f < F < 1.8, TcrSss grows as a power-law with the size L, as shown here in log-log scale (see Eq. [30}. (b) Corresponding 
anomalous exponent n(F) defined in Eq. [30] as a function of the force F. 
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FIG. 5: Critical exponents in the limit of vanishing force F — ¥ (the points correspond to F = 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4) 
(a) Log-log plot of the exponent fi(F) as a function of F : the slope corresponds to a ~ 0.6 (see Eq. I32|l (b) Log-log plot of the 
correlation length £(F) as a function of F : the slope corresponds to v ~ —1.29 (see Eq. I33[) 
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As shown Fig. [4] (a) , we find the following typical behavior in L 

1 



lriT crossmg (F, L) ~ (lni - ln£(F)) (30) 



where the corresponding exponent n(F) is shown as a function of the external force F on Fig. [4] (b). 

After inversion, the behavior of the typical crossing time T cross i ng (F, L) corresponds to the following dynamical 
behavior for the typical displacement x{t) in the direction x of the force, as a function of time t 

Xtyp(t) ~ t(F) (31) 

Our results are thus compatible with a zero- velocity phase characterized by an exponent < fi(F) < 1 in a finite 
region < F < F c , and with a finite velocity phase with n(F) = 1 above the threshold F > F c . We cannot give a 
very precise value for the threshold F c as a consequence of the rounding coming from finite-size effects. Anyway, this 
critical value is not universal and depends on the details of the random potential at short distances. As a consequence 
in the following, we will focus instead on the universal properties that are expected to appear at small external force, 
where the dynamics is governed by the large-scale properties of the random potential. 

In the region 0.1 < F < 0.4, our results are compatible with the power-law 

u(F) oc F a with a ~ 0.6 (32) 

F->0 

as shown on Fig. [5] (a). It is thus significantly different from the exactly known value a(d = 1) = 1 for the same 
problem in dimension d — 1 (see the review [J] and references therein). Similarly, as shown on Fig. [5] (b), we find 
that the correlation length £(F) appearing in the prefactor of Eq. [22 (that we measure from Eq. 1301 by fitting our 
data shown on Fig. @] (a) diverges with the following power-law as F — > 

f(F) oc F~ v with v ~ 1.29 (33) 

Again this exponent is thus significantly smaller than the exactly known value v(d = 1) = 2 for the same problem in 
dimension d — 1 (see the review Q and references therein). 



V. DISCUSSION : EQUIVALENCE WITH AN APPROPRIATE RENORMALIZED DIRECTED TRAP 

MODEL 

A. Scaling analysis of the equivalent problem in dimension d — 1 

As recalled in the introduction, the equivalent problem in dimension d = 1 is exactly soluble by various methods. 
However, since these exact methods cannot be applied in dimension d = 2, it is important to identify the physical 
mechanism of the anomalous diffusion in dimension d — 1 to see if scaling arguments can be used in higher dimension 
d > 1. In dimension d = 1, the physical picture is that the dynamics renormalizes onto a directed trap model [4l.l9l.IT0j 
with the following properties : 

(i) the distribution P(t) of trapping times decays as the power-law 

P ^K^^m (34) 

where the exponent n(F) is exactly the exponent governing the anomalous diffusion law of Eq. 1311 Equivalcntly, the 
distribution P(B) of barriers B = lnr against the drift F decays exponentially 

V(B) oc e -^ B (35) 

This means that fi(F) should be interpreted here as the inverse of the typical barrier 

B t y P (F) = J- (36) 
^ ) 

(ii) the size £(F) of a renormalized trap scales as 



(37) 
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as a consequence of the scaling property of the random potential in terms of the Hurst exponent H that yields 
Btyp(F) = (£,(F)) H . So the two exponents governing the divergences of barriers 



and lengths 



£(F) oc (39) 



are related by the relation 

v = | (40) 

(iii) the exponents a and v governing the divergences of the typical barrier and of the typical size in the limit of 
small external drift F can be understood via the following simple scaling argument in dimension d = 1 : the barrier 
B(l) against the drift F on a interval / can be written as 

B(l) = -Fl + l H u (41) 

where u is a random variable of order 0(1). The minimization of B(l) with respect to I yields the characteristic length 
scale 

^ ld \F) cx ^ with uu = (42) 

above which the deterministic term {—Fl) becomes larger than the random term in Eq. 1411 So £} ld \F) represents 
the length of an effective renormalized trap, and the corresponding barrier scales as 



R (W) 



(F) = B(l = e d) (F)) cx (e d) {F)) H oc JL W ith a u = ^ = Hv ld (43) 



In dimension d = 1, the precise equivalence between the biased Sinai model and the one-dimensional directed trap 
model is discussed in more details in In particular, besides the scaling of Eq. [3_TJ one can relate exactly the 

corresponding diffusion fronts that involve Levy stable distributions. We refer to [10| and references therein for more 
details. 



B. Similarities and differences in dimension d = 2 



In dimension d = 2, the numerical results presented above suggest that, as in dimension d = 1, there exists an 
anomalous diffusion phase as small external force, which is equivalent with a directed trap model at large scales. The 
important difference is then in the properties of the elementary effective renormalized traps, in particular in their 
dependence upon the external force. As explained in (i) and (ii) above, fi(F) should be interpreted as the inverse 
of the typical barrier of an elementary trap, whereas £(F) should be interpreted as its longitudinal length in the 
direction of the force. Taking into account the numerical uncertainties on the the exponents a and v measured in 
section HVl the values given in Eqs I3"2l and 1331 seem still compatible with the scaling relation of Eq. [37] with H = 1/2. 
Our conclusion is thus that only the point (iii) which contains an explicit one-dimensional optimization is not valid 
anymore : we are not aware of the appropriate argument based on a two-dimensional optimization that would give 
values of a and v in d — 2. The fact that these traps are 'smaller' as a function of F than in the equivalent problem in 
dimension d = 1 means that the particle uses the transverse direction to find lower barriers. However, a quantitative 
analysis of this effect goes beyond the present work. 



VI. CONCLUSION 



In this paper, we have studied the continuous-time random walk of a particle in a two-dimensional self-affine random 
potential of Hurst exponent H = 1/2 in the presence of an external force F. We have presented numerical results 
on the statistics of escape probabilities and of first-passage times that satisfy closed backward master equations. Our 
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main conclusion is that, as in dimension d = 1, there exists a zero- velocity phase in a finite region of the external force 
< F < F c , where the dynamics follows an anomalous diffusion law x(t) ~ £,(F) t^ F ' with < fi(F) < 1. However, in 
the limit of small external force F — > 0, the exponents governing the vanishing of the anomalous exponent n(F) oc F a 
with a ~ 0.6 and the divergence of the correlation length oc F~ v with v ~ 1.29 are significantly different from 
their exactly known values a(d = 1) = 1 and v(d — 1) = 2 in dimension d = 1. Finally, we have proposed the following 
interpretation : as in dimension d = 1, the dynamics renormalizes at large scales onto a directed trap model, so that 
the length and the inverse of the anomalous exponent represent respectively the typical length along 

the direction of the force, and the typical barrier l/fj,(F) of an elementary renormalized trap. The fact that these 
elementary renormalized traps are 'smaller' in linear size and in depth than in the equivalent problem in dimension 
d = 1, indicates that the particle uses the transverse direction to find lower barriers. An important issue for the future 
is thus to analyse more quantitatively this effect in order to better understand the values of the exponents a and v. 

Another way to understand these values could be via a strong disorder renormalization analysis : we have shown 
in our previous work [i~3j that without external force F = 0, the dynamics is governed by an 'Infinite Disorder 
Fixed Point', so we expect that in the limit of small external force F — > 0, the dynamics should be described by a 
'Strong Disorder Fixed Point' as in dimension d — 1 [11]. It is thus interesting to make some comparison with the 
two-dimensional Random Transverse Field Ising Model (RTFIM) which has been much studied via strong disorder 
renormalization [27l - l30j (see also the Monte-Carlo study [3l| on the RTFIM and the Monte-Carlo study [32| on the 
two-dimensional disordered contact process that have found evidence for infinite-disorder behaviors) : it turns out 
that the most recent values of the correlation length exponent are around vrtfim — 1-25 [29l I30L l32j. which is rather 
close to the value v ~ 1.29 that we have measured here. Although the RTFIM and the Sinai model are completely 
equivalent in dimension d = 1 [TTj . the two-dimensional RTFIM and the dynamical model considered in the present 
paper are not directly related : in particular, in the RTFIM, the disorder is local, and the problem remains an isotropic 
two-dimensional model, whereas in the dynamical model, the random potential presents long-ranged correlations, and 
the external force introduces a very strong anisotropy since the model becomes essentially directed in the direction 
of the force. From this point of view, the fact that similar values for the exponent v have been found seem to be 
purely a coincidence. On the other hand, it turns out that the value H = 1/2 of the Hurst exponent of the random 
potential that we have chosen to define the dynamical model, is the analog of the exponent if) ~ 0.5 measured in 
the two-dimensional RTFIM (2|| H3, [HJ ■ So a more detailed comparison between the strong disorder fixed points 
associated to both types of models is needed to determine whether the corresponding exponents are close just by 
coincidence or not. 
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